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Abstract 

We derive an expression for the collective Casimir- Polder interaction of a trapped gas of condensed 
bosons with a plane surface through the coupling of the condensate atoms with the electromagnetic field. 
A systematic perturbation theory is developed based on a diagrammatic expansion of the electromagnetic 
self-energy. In the leading order, the result for the interaction-energy is proportional to the number of 
atoms in the condensate mode. At this order, atom-atom interactions and recoil effects lead to corrections 
compared to the single-atom theory, through shifts of the atomic transition energies. We also discuss the 
impact of the spatial delocalization of the condensate mode. 



1 Introduction 

It is well known from cavity quantum electrodynamics (cavity QED) that the energy levels and lifetimes of 
the electronic states of an atom placed near a macroscopic body are shifted from their free-space values l|2 . 
This effect can be understood from the modification the body imposes on the vacuum field modes which 
lead, for example, to a position-dependent change in the Lamb shift. The resulting (van der Waals or 
Casimir-Polder) force between the atom and the macroscopic body has been shown to match the predictions 
of QED in several experiments [^|^ ■ 

As cavity QED effects often do not require a relativistic treatment of the electronic or atomic motion, 
the techniques traditionally employed are lent from non-relativistic QED: a mode expansion of the elec- 
tromagnetic field and a first quantized theory for the remaining (atomic) part of the system [8j. Instead 
of working with mode expansions adapted to the presence of a body, there is another approach making 
use of the fluctuation-dissipation theorem (9j[l0]: the level shift is cast in a form involving generalized 
susceptibilities from linear response theory, the (retarded) Green functions. The influence of the surface is 
then encoded in the appropriate scattering amplitudes of the body, e.g., reflection coefficients for a planar 
interface. This makes the approach applicable to very general descriptions of the surface material, including 
absorption and dispersion. Another advantage of the formalism lies in the fact that renormalization gets 
simplifled, as the (divergent) free-space part of the Lamb shift is easily isolated from the surface-dependent 
contributions, the latter being finite. 

In the present paper, we are interested in the shift of the collective energy levels of an A^'-atom system 
due to the presence of a nearby surface. A theory that has to account for the quantum statistical character 
of atoms is conveniently formulated in terms of second-quantized atom field operators. We follow the 
standard procedure for perturbation theory, which offers a pictorial representation in terms of Feynman 
diagrams and permits us to calculate the elements of the electromagnetic self-energy, approximating them 



with the Dyson series |llj. For our purposes, the theory has to deal with a confined atomic system in a 
trap (including inter-atomic interactions) , and the interaction with the electromagnetic field is the relevant 
perturbation. While much of the literature on Bose-Einstein condensates (BECs) in an external potential 
deals mainly with the collective properties of atoms in their electronic ground state, a quantum field 
theory of ultra-cold atoms interacting with photons was formulated in 12 13]. We build on this approach 
and merge it with the linear response techniques for electromagnetic field fluctuations near a surface. 
In the present paper, we consider the atom-light interaction up to second order, which is the flrst non- 
vanishing contribution. We find under quite general circumstances the atom-surface interaction energy and 
demonstrate that it does not reduce to an integral over the density distribution of trapped (ground state) 
atoms. The propagation in the excited state, although only virtually, connects ground state correlation 
functions at different space-time points. This leads to a recoil shift of the atomic polarizability, in addition 
to the familiar density shifts due to the atom-atom interaction. For the atomic ensemble, we consider two 
simple examples: firstly, an interacting BEC at temperatures well below the critical temperature, where 
we consider only a single mode of the atomic field with macroscopic occupation. Our second example is 
the ideal Bose gas at nonzero temperature that can be essentially characterized analytically. Both systems 
are held in harmonic traps centered near the surface. We develop in this paper the main methods, check 
that several limiting cases are recovered and discuss the two examples above in some detail. The aim for 
future publications is to generalize this approach in two respects: on the cavity-QED side of the problem, 
to push the atom-field interaction to higher orders and, on the BEC side, to take into account low-lying 
collective states of the interacting atomic ensemble like Bogoliubov quasi-particles. 

The paper is organized as follows: In section [2] we describe the interaction of the atomic system and 
the electromagnetic field, in a form involving second quantized operators for atoms as well as for the 
field. In section [S] we calculate a general expression for the energy shift of the atomic ensemble due to 
the interaction with the electromagnetic (e.m.) field by evaluating the first non-vanishing term in the 
Dyson expansion of the S-matrix. The result thus obtained involves the Feynman propagator for the e.m. 
field in the presence of a surface which is introduced in section |4.1[ Atomic propagators are calculated 
in section 



4.2.1 



and section 4.2.2 for an interacting BEC (zero temperature) and an ideal gas (nonzero 
temperature), respectively. The results of section [2] to section [4] are then used to calculate the atom-surface 
interaction of these two examples (section |5] and section [6|. We cross-check our calculations against existing 
results in section [sTil by re-deriving the Casimir-Polder potential for a single perfectly localized atom. 

Our units are such that h = ks ~ 1, the speed of light c and the atomic mass M are kept for the ease 
of reading. 



2 Quantum field theory of atoms and photons 

We consider A'' identical atoms in a trap above a fiat surface. The surface is taken to lie in the xy-plane, the 
center of the trap is located a distance d from the surface in the half-space 2: > 0. The atoms are treated 
in the electric dipole approximation with an electric ground state \g} and excited states |e). The extension 
of this model to more realistic atoms is straightforward by summing the contributions of all excited states 
in the calculation of the ground state shift. 

Apart from possible inter-atomic interactions, the atoms interact with the electromagnetic field via a 
fj, ■ E interaction term, where the dipole operator has transition matrix elements /x®*^ — {g\d\e). (For a 
comparison between the minimal coupling Hamiltonian and /x- E interaction, see |14[|15| .) The interaction 
between the atomic system and the surface originates in this atom-field coupling: the surface contains 
sources that radiate a field, and it imposes boundary conditions on both the intrinsic field fluctuations and 
the fleld radiated by the atom. The relevant correlation functions of E near the surface will be dealt with 
in section |4T] 

As in |12|, we will work with a Hamiltonian that describes the atomic degrees of freedom (as well as the 



2 



electric field) in second quantization i.e., a quantum field theory of atoms interacting with photons. The 
operators *I'g(r) and ^'e(r) describe the annihilation of an excited-state or ground state atom at location 
r. 

As we want to treat the influence of the electromagnetic coupling as a perturbation to the atomic 
system, we split the total Hamiltonian as follows: 

H = Ha + Haf+Hf (2.1) 

Here, Hp is the Hamiltonian for the unperturbed field in the presence of the surface, Haf contains the 
atom-field interaction, and the Hamiltonian Ha describes the trapped atoms. The atomic operators in an 
interaction-picture with respect to Haf then have the general form 

= (2-2) 

n 

*^(^) = / exp(iq.r)e„(t) . (2.3) 



where the time dependence of the operators g-n{t) and e^it) is specified in section 4.2 below. They satisfy 



the bosonic or fermionic equal-time commutation relations. In terms of these field operators, the atom-field 



interaction Haf in eqn. (2.11 can be written as 



Haf = - j d^a;^{i5.(x)[Mr*I(a;)*e(x) + M:''*i(a;)*sW]} (2-4) 



(compare 12 eqn. (81)]). As mentioned above, it is this term that is responsible for the interaction between 
the surface and the atoms, as the specific form of E(x) depends on the surface. We do not make the rotating 
wave (or resonance) approximation here because otherwise relevant virtual processes would be missed. 

We use the notation r = (x, z) for spatial vectors, where the two-dimensional vector x lies in the plane 
perpendicular to the surface. Spatial integrations J (f'r run only over the z > Q half-space. Spacetime 
points are denoted by a; = (r,t). 



3 Second-order energy shift 



The aim in this section is to calculate the energy shift of the atomic system due to its interaction with 
the electric field. In the case of a single atom in front of a surface, this shift is usually calculated in 
time-independent perturbation theory [9l |10[[l6] . We will employ instead standard tools from field theory: 
the energy shift is obtained from the S-matrix, which can be perturbatively approximated with the Dyson 



series (see 17 sec. 3.5]). For a treatment of the single atom in front of a surface in both formalisms, 
nonrelativistic perturbation theory and the Dyson series, see [18||19| . 

Let us briefly recall the basic relations which will be used: the energy shift of an unperturbed state of 
the atomic system can be calculated from the real part of the self-energy (logarithm of the S-matrix). In 
the present paper, we will consider only terms up to the second order in Haf in which the self-energy and 
the T-matrix coincide. Recall that the T-matrix is deflned as the nontrivial part of the S'-matrix, 

Ss^ = 5{f - i) - 2-Ki S{Ef - Ei) Tf, , (3.1) 

which, in turn, can be expressed as a series of time-ordered products of interaction picture operators, the 
Dyson series: 

_5^^_|_^(z*r f dtl...dtnT{HAF{tl)...HAF{tn)} , (3.2) 

where the symbol T{. . . } denotes time ordering. 
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For a general self-interacting atomic system, it is convenient to define the interaction-picture operator 



*(2;) = ft (a;) 4-^(3;) 



(3.3) 



as the operators and 'J'e appear only in this combination or its hermitian conjugate in eqn. (2.41. With 



the initial and final states containing no excited-state atoms, the first-order term in eqn. ( |3.2[ ) vanishes, 
leaving the second-order contribution 



5(2) 



eg qe , 



(3.4) 



(3.5) 



The brackets (...) in eqn. ( |3.4[ ) denote an expectation value in a stationary state of the atomic Hamiltonian 
Ha- In the above diagram, the in- and outgoing lines represent A'^ atoms in the state \g) that make up the 
unperturbed atomic state. The virtual state (inner line) consists of an atom in the state |e) (dashed line) 
propagating in the presence of a background field (solid line) made up of the remaining A'^ — f ground state 
atoms (still a large number). The vertices, where an excited atom is created or destroyed, are proportional 
to the dipole moment of the transition: 



i -en ^ ■ <7e 
& - - — —lHa , - - O ~ —l-^-a ■ 

The photon line in eqn. (|3.5|l is given by the time-ordered (or Feynman) propagator 



Dal3{xi;X2) 



{T{E^(x-,)Efi{x2)}) 



duj 
2^ ' 



(3.6) 

(3.7) 
(3.8) 



where the brackets (...) in eqn. (3.71 denote an expectation value with respect to an equilibrium state of 
Hf. 



Finally, the contraction in eqn. (|3.4[| is defined as 



^{X2)'f\xi) = [*(X2), *^(a;i)] e(t2 - h) . 
which can be decomposed for bosonic or fermionic fields (upper/lower sign) as 



(3.9) 



*(x2)*t(a:i) = e{t2 - h){±[¥,{x2),^Uxi)]^ [*e(a;2), 

+ <5flix2) [*e(a;2), ■^lixi)]:f ^ gixi) ± *i(xi) (^2), •fs{xi)]^ e{x2) 

+ ['i'l{x2),'i'l{xi)]:f<lig{xi)<l'4x2) + *3:(a;i)*^(a;2)[*e(a;2),*9(a;i)]T} 



(3.f0) 



If our initial and final states contain no excited atoms, the last three terms will yield zero in an expectation 
value, and we are left with 



{<f{x2)<S'Hxi)) = (*J(x2)*e(2;2)*i(xi)*9(xi))e(t2 " h) = 



(3.11) 



for both statistics. 



We will see in eqn. (4.171 below that for an ideal gas, the above expression reduces to the form that 



is usually obtained from applying Wick's theorem to a time-ordered product of four interaction picture 
operators (see [11[ chap. 3]). This is no longer true in the general case (interacting atoms), and the 



Feynman-rules for translating a diagram containing a line like eqn. (3.111 must take into account the 
presence of other lines due to the interaction with the background field. 
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4 Photon and atom propagators 



function Df^g of eqn. (|3.7l 



In order to evaluate the general expression eqn. ( |3.4[ ), we now need to assume a concrete form for the 

the propagator of the electric field in the presence of a surface — and for 



the expression (^g(a;2)^'e(a;2)*I'e(a::i)*g(a;i)) in eqn. (3.111, characterizing the atomic ensemble. For the 



former, we can rely largely on work presented in 110 , which will allow us to apply our technique to very 



general surface materials. Concerning the latter, we will focus on a pure condensate in a trap (section 4.2.1 1 
and on a trapped ideal Bose gas at nonzero temperature (section |4.2.2[ |. 



4.1 Photon propagator near a surface 



The time-ordered propagator for the E-field [eqn. |3.7l] is usually worked out explicitly from a mode 
expansion of the E-field. This can be done in the presence of a non-dispersive surface, too, with the mode 
functions getting of course more cumbersome to satisfy the boundary conditions at the surface |19[|20| . 
We want to follow here the approach of [9l |10[[2l] , which connects the field propagator to a form involving 
correlation functions from linear response theory, the retarded Green functions. This applies as long as the 
fluctuation-dissipation theorem for the electromagnetic field holds [21[|22| . 
The retarded Green function for the electric field is defined as 



Gal3ixi,X2) 



(4.1) 



with Fourier transform Gai3{ri,r2,oj). By rearranging the time ordered product in eqn. (3.7 1 and using 



the fiuctuation-dissipation theorem (see |10( appendix B]) we can express the Fourier transform of the 



Feyiunan-propagator eqn. (3.81 as 



Df^(ri,r2,aj) = Im [GQ^(ri, r2, a;)] coth[ — ] - iRe [Ga0(ri,r2,ui) 



(4.2) 



(see [11[ sec. 31]), where T is the temperature (fcs = 1) of the field. Here we assume the field and its 
sources in thermal equilibrium at the temperature Tp. The atomic part of the system may have a different 
temperature and is even allowed to be in a non-thermal state. As we will see below, it is preferable to 
integrate the retarded Green function along the imaginary frequency axis. Using the fact that Ga0{oo) has 
only poles in the lower half of the complex tj-plane, we can express the cj-integration in {No\T^^^No) (see 



eqn. (5.11 below) as (see 16 Appendix A]) 



/ 



duj 



oJq — — i 



d^Gc.fi{iO 



27rGQfj(— cJo)6(— cjo) 



(4.3) 



where cjq 



,(q, TVo). We have made here the approximation iuo 2> Tp, i.e., field temperatures much 



smaller than the atomic transition energies, where the summation over the poles of coth([j/2r) can be 



replaced by an integral. The second term on the rhs of eqn. (4.3 1 is nonzero only for excited state atoms 
{Ujq < 0) and describes spontaneous emission and resonant contributions to the energy shift 16 23 . For 



atoms in the ground state, corrections to Eq.(4.3 1 are proportional to the number of thermal photons which 



is exponentially small if Tf is much smaller than the relevant transition frequencies. When using eqn. ( 4.3 1 



in the remaining sections, we will suppose throughout that the number of thermal photons is negligible, 
and any temperature dependence that appears from now on is always associated with the temperature of 
the atoms, not the photon field. The generalization to finite field temperatures is left for future work. 

Now, from linear response theory (see 11 sec. 32]) and the linearity of the Maxwell equations, the 
response function Gc«^(ri, r2, cj) can be identified with the classical Green function, i.e., the electric field 
at ri generated by a classical dipole, oscillating at frequency lj, which is located at r2. The explicit form 



of the Green function in the presence of an interface is well known 24 and can be split into a free space 
and a reflected part: 



GaB — Gf^d + G, 



'all 



(4.4) 
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where 6*% is the retarded Green function in free space. As we are only interested in that part of the energy 



shift caused by the presence of the surface, we will not consider at all. The decomposition eqn. (4.4 1 
permits us in a simple manner to subtract the divergent diagrams involving photon loops that yield the 
free-space Lamb shift, because the latter arise from the Green function G°^. To get the distance-dependent 
part of the energy shift, we will simply substitute Gai3 by G^^. The expressions containing G^^ are then 
finite without any further renormalization. 

The surface contribution G^^ at imaginary frequencies has the form 

G«^(ri,r2,iC) = / — i?.;3(C,k)e-«(^i+^^)e"-(-i--^' , (4.5) 

ItT J K 

(see [9||24| and below in appendix |a] for more details) where fio = (eoc^)~^ is the vacuum permeability and 
K = \/(^ jc? -f fc-^. The two-dimensional vectors x and k denote position and momentum vectors parallel 
to the surface, respectively. The tensor elements R^p contain the reflection coefficients appropriate for the 
specific surface material. As we are only considering the reflected part G^^, we will in the following skip 



the label 'R' from eqn. 1 4.5 1. Note that from the viewpoint of perturbation theory, the surface response 
functions R^p depend on the quantum state of matter in the surface; they are calculated, of course, in the 
absence of the atomic system outside it. 

4.2 Propagators for atoms 



In the following, we calculate the atomic part of expression eqn. (3.41 for two simple examples of atomic 
systems. Together with the photon propagator obtained in section [4T] above, these will flnally allow us to 
evaluate the atom-surface interaction in section [5] and section |6l 

4.2.1 Dilute interacting BEC in the single mode approximation 

For the interacting dilute Bose gas confined in a trap, we further restrict ourselves to the deeply degenerate 
case, where we can consider a large number A'^o of atoms in a single condensate mode. We leave the 
contribution of condensate (Bogoliubov) excitations for future work. The atomic Hamiltonian Ha describes 
two-level atoms with a contact interaction between excited and ground state atoms: 



Ha 



E(gogo) + J d^q {'^lini'i) + bgeglgo) e^eci (4.6) 



Here, the energy a;e(q) — uje + q'^/2M contains both the electronic excitation energy and the kinetic 
energy. The constant bge characterizes the interaction between ground- and excited state atoms. The 
self-interaction amongst the ground state atoms and the effects of the trapping-potential are contained in 



E^glgo). In our approximation, the field operator 'ifg in eqn. (2.2 1 consists only of a single mode with 
the mode function <^o(r), which is the condensate wavefunction calculated self-consistently by solving the 
Gross-Pitaevskii-equation (GPE) [25] 

Here, j-i{N) = dE{N)/dN denotes the chemical potential, the constant bgg characterizes the self-interaction 
of ground-state atoms, Vtg denotes the trapping potential felt by the ground-state atoms, and the condens- 
ate wavefunction is normalized to / (fir |0o(r)P — 1. 

Interactions between excited state atoms are neglected in this paper. This is legitimate since our 
unperturbed state consists of a large number of ground state atoms. Excited state atoms will then only 
occur in virtual states, and their number will be small. 
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With the particular choice eqn. (4.6 1 for Ha, the field operators 'i/g and from eqn. ( 2.2 I and eqn. ( 2.3 1 
assume the form 



= (l)o{r) exp[-'i{Eg{No + l)-Eg{No) + bgeNe)t]go , 



(27r)3/: 



exp[i (q.r - (tJfein(q) + bge No) t)] 



(4.8) 
(4.9) 



The diagram eqn. (3.111 then becomes 

^—^^ = e{t2 - ti) <j}o{ri)<t>o{r2) No 



^3 

(27r)3 



-il^cg(ci,No)(t2-ti) 



Here the transition frequency Ueg (q, A'^) is defined as 

UJ,g{(l, N) = a;e(q) + Eg{N - 1) - £g(iV) + (iV - l)6ge , 



(4.10) 



(4.11) 



where the frequency shift of the atomic transition due to inter-atomic interactions appears. If the system 
consists of a single atom only, 

„2 



a;eg(q, 1) = cje(q) 



(4.12) 



- "^^^ ■ 2M 

which is the resonance frequency of a single atom, including the recoil shift. The physical interpretation of 
Eq. (4.101 is quite clear: a virtual photon takes a ground state atom at position ri to the excited state, the 
atom propagates freely to position r2 and joins the other ground state atoms there. We shall see below that 
the relevant distances \v2 — ri| are negligibly small so that eventually the ground-state density |0o(ri)P 
determines the atom-surface interaction. 



4.2.2 Ideal Bose gas at finite temperature 

For the non-interacting trapped Bose gas, treated in the grand-canonical ensemble with a mean total 
particle number N , an inverse temperature P and chemical potential /x, the Hamiltonian Ha takes the 
form 

Ha = Y.^^ + I d'9^Ln(q)4eq • (4-13) 
The mode functions for the operator ^ g are the single-particle wavefunctions 0n that solve 



h^ + ^t.W] <^n(x) = iJ„<^„(r). 



The field operators eqn. (2.2 1 and eqn. ( |2.3[ ) now take the simple form 



(r)e " Qn 



(27r)3/2 



The atomic part eqn. (3.111 yields 

Af, /!, /3 ■ 'N. iJ.il = 



(*t(x2)*,(Ti))vI/,(a;2)*l(a;i) 



(4.14) 

(4.15) 
(4.16) 



(4.17) 
(4.18) 



X y (^jj^s exp[iq.(r2-ri)-i((7V(2'n)-fcJe9)(i2-ti)] 

Note again the occurrence of the two-point correlation function for the ground-state atoms. We thus reach 
a similar structure as in Eq.(4.10l above, but with a sum over all trap eigenstates. 
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5 Energy shift of an interacting Bose gas trapped near a 
surface 



5.1 Generalized polarizability 

With the results obtained above, we can now evaluate the interaction potential between the single-mode 



condensate and a surface. Putting the expression for the photon propagator eqn. (3.7 1 and for the atomic 
two-point function eqn. (4.101 into eqn. (3.4 1, we get for the T-matrix element (after performing the dti 
and dt2 integrations) 



(iVo|r^''|7Vo) = No/ii'fil' I dVi / dV2 0o(r2)<?!'S(ri) 



_^ o»q'(ri-r2) 

X / T^-DQ^(ri,r2,a;) , r ^r ^ ^ • ' 

271*^ J (2tv)'^ oj - ujcg(ci, No) + le 



(5.1) 



Using relation eqn. ( 4.3 1 in eqn. (5.11 to link the time-ordered photon propagator to the Green tensor Gai3 , 
we get 

2 



■•2) 



{No\T^^^\No) = -^oj^^i't^;' J rf'ri J dV^ <^o(r2)<^S(ri) 
X d^G.;3(ri,r2,iC) J d^ga(q,e)e''> (--i 

The generalized polarizability 

^e,(q,jVo) 
"^'^'^^ <(q,7Vo)-fC^ ■ 
contains the interaction- and recoil-shifted resonance frequency uugg{q, Nq) (see eqn. (4. 111). In eqn. (5.2l 
we neglected the resonant contribution of thermally excited photons. 



(5.2) 
(5.3) 



5.2 Condensate wave function 



In order to evaluate eqn. (5.2 I, we have to substitute a suitable approximation for the condensate wave 
function (?!)o(r). For simplicity we solve the Gross-Pitaevskii equation eqn. (4.71 with an isotropic harmonic 
trapping potential 



■ (x' + {z- df) 



(5.4) 



where d denotes the distance of the trap center from the surface. If the kinetic term in the GPE can be 
neglected (Thomas-Fermi approximation), the solution for the density profile takes the form of an inverted 
parabola. This is usually a good approximation for large particle numbers. Here, we choose a Gaussian 
ansatz for the wave function because it simplifies the subsequent integrations. (For calculations with a 
Thomas-Fermi profile, see 26 .) The ansatz also allows for the limit TVq — >■ 1 in order to provide a cross- 



check with results for a single-atom system (section 5.41. Gaussian functions also approximately solve the 
GPE, if width and amplitude are varied such that the Gross-Pitaevskii functional is minimized (see 25 
for details). We thus make the ansatz 

' 2aHNo) \ ■ ^^-^^ 
The minimization procedure gives a spatial width a{No) in eqn. (|5.5|l that depends on the number of 



(^o(r) = {y^a{No)) ^ exp 



trapped particles and has the asymptotic values 27 



aiNo) 



oo, 
ao 



No = 1 
2"iVoay''^ iVoffl 
■K ao ) ' ao 



> 1 



(5.6) 
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where ao = {M v)-^l^ is the width of the single-particle ground state in the trap. The s-wave scattering 
length a is related to the interaction constant hgg from eqn. (4.7 1 via 633 = Ana/M. In the second case of 
eqn. ( |5.6[ l, the interaction energy of ground state atoms is much larger than the bare harmonic potential. 
This regime corresponds to the Thomas-Fermi limit (the Thomas-Fermi radius is Rtv ~ ao(15Afoffl/ao)^''^). 

A subtlety arises for the Gaussian ansatz (5.5 1 because it is normalized only in the limit d ^ cr(A^o) 
if spatial integrations are restricted over the half-space z > Q. We shall always assume this limit, as 
our approach is clearly not valid for atoms touching the surface. The wave function <^o is of the order 
0{ey^p[—{d/ (j(Nq)Y\) at the surface, and exponentially small terms of this order will be systematically 
discarded in numerical evaluations of energy shifts in section |5.4| and section [6] These approximations are 
dealt with in detail in appendix [B] 

In the following, we will evaluate ^jvoVq with the approximate ground-state wave function 00 from 
eqn. (5.51. The result eqn. (5.2 1 for the T-matrix is, however, more generally valid and can be evaluated 
similarly for other approximations of (/>o(r). 



5.3 Recoil shift and (de) localization correction 

At this stage, it is convenient to introduce sum and difference coordinates r± and to split them in com- 
ponents perpendicular and parallel to the surface: z± = Zi ± 22, and x-t = Xi ± X2, Similarly, for the 
momentum, we use from now on q = {qx,qy,qz) — > (q, 9z)- Integrating in eqn. (5.21 over ri and r2 and the 
angle of the two-dimensional vectors k and q (k • k = fc^ and q • q = g^), we get 



(ATolr'^'liVo) = 



T-3/2 



eo 



K 2 



qdqa{q,q^,^) Io[2cr'kq]t 



- e" - (1 -I- erf [ na]) Map 

a 



(5.7) 



with an obvious notation for a(g, Qj,,^). The diagonal matrix Map originates from the scattering tensor 
Rap and has elements 



M,4fc,e) = Myy{k,i) = R''k^ + {R" - R') {i/cf 
M,,(fc,e) = 2R''k'^ , 



(5.8) 
(5.9) 

and the i?^'^(fe,^) are the reflection amplitudes from the surface (eqns. ( |A.3[ | and ( |A.4[ )), Jo is a Bessel 
function of the second kind, k — sjk"^ + and erf denotes the error function. 

To perform the dq^, and dg-integrations in eqn. (5.71, we observe that in a(q,qz,i,) (see eqn. (5.3l) the 
momenta q and appear only as recoil shifts of the atomic transition frequency uieg{No) (see eqn. ( 4.11 1). 
Since the relevant momenta are limited to typically 1/cr, the recoil shift is a small correction because 
l/(Mo-^) = V <^ ujeg{No) is usually well satisfied. We therefore expand in powers of q and q^ and integrate 
term by term by means of the identities 



Jo 



qdqe-^'' /o[2fcga" 
qdqq e ^ ^ ' Iq [2kqa 



1 

2^ ' 
— (- 



Thus, we finally obtain for the T-matrix eqn. (5.7 1 



(ATolr'^'lAo) = - 



Ao fJ-i'^fJ-f 
2n £0 



kdk 1 

^fiT 2 



+ k' 



'(l + erf[- -Ka])Mai3 
a 



(5.10) 
(5.11) 



x{a{^,No) + a^'-'\^,No,k)} 



(5.12) 



9 



where the polarizability 



describes the no-recoil case. The recoil term a''^'^' is given by 



(rc) 



+ 



^4M(cr(iVo))2 ' 2M J 
We can attribute this correction to a recoil shift of the effective resonance frequency 



+ 



4:M{a{No)y 2M 



(5.13) 



(5.14) 



(5.15) 



where the two terms describe the kinetic energy from the delocalized condensate wave function and from 
the absorbed photon momentum in the exc ited s tate, respectively. 

(2) J l 

The T-matrix element Tj^^j^^ from eqn. ( 5.12 1 is our main result for the interaction energy of a trapped 

f 21 

Bose gas with a plane surface. In the above form, it is clear that Tj^^j^^ generalizes the result for a stationary 
single atom in a straightforward manner. Clearly, as we put A^o = 1, we get the single-atom transition 
frequency tJeg(l) = cueg- And with the identity 



hm - e 

cy-fO 2 



(1 -h erf [- - H) = e" 



we get from the no-recoil term of eqn. | |5.12[ ) 



hm r(f 



qe eq 



+e 



(5.16) 



(5.17) 



with ro — (0, 0, d) the position of the trap center. This is the known result for a perfectly localized 
single atom as in Ref. 10. eqn. (2.28)]. The recoil correction involving a'""' is discussed in more detail in 
section |5.4[ It is usually very small, unless the trap frequency u is comparable to the atomic resonance 



LUeg, a case of no practical significance. 

For a large atom number No, the resonance frequency ijJeg{No) in eqn. (5.131 incorporates the inter- 
atomic interactions (see eqn. (4.11 1). The overall proportionality factor A^'o of eqn. ( 5.12 1 can be understood 
by recalling that the responsible diagram (see eqn. (3.5 1) represents a sum of self-energies of A'^o individual 
ground state atoms. In higher orders, i.e., diagrams with four or more vertexes, virtual photons can connect 
different ground state atoms, and we can expect a nonlinear scaling in A'^o- 



5.4 Single ground-state atom 



The T-matrix element for a single atom can be obtained from eqn. (5.121 by setting A^o ~ 1- Introducing 
the scaled distance x = dueg/c, and rescaling the integration variables ^ — ^/ujeg, k = ck/ueg, k = cn/ujeg, 
the T-matrix reads 

(liT'^'lD = ^^f^ Td? r^/(7.,.,r;)M„,(fc,?) 

XLJeg {a!(Jajea, 1) -f a'"''(Ci^es, 1, fcoJej)} , (5.18) 

where the energy scale is set by the natural linewidth jeg ~ |/i*'^pajfg/37reoc"^, the so-called Lamb-Dicke 
parameter rj = u^gao/c gives the size of the trap ground state in units of the resonant wavelength. The 
quantity / becomes 

1 X 

I{k, x,ri) = - exp[— 2kx -|- «: ] (1 + erf [ Tt?;]) (5.19) 

2 T) 
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The matrix Ma/s defined in eqn. (5.91 depends on the reflection coefficients and R" and encodes the 
surface properties. In the dimensionless units of eqn. (5.181, the recoil correction a^^"^^ is now seen to be 
proportional to the ratio f/cjeg'- 



UJeg a{^ijJeg, 1) 
X^"'"^ (^i^cg, 1, kljJe.g) 



V 1 — C /'3 k 2\ 



(5.20) 
(5.21) 



The trapping frequency v/2'k for a single ground state atom in the potential eqn. (5.4 1 is usually around 
10 . . . 1000 Hz, much smaller than the frequencies of optical transitions ~ IQ^^ Hz. This justifies the 

expansion of the recoil shift for small atom momenta q done in section[5] Experimental situations where the 
recoil correction is enhanced in magnitude could involve tight traps like optical lattices {p/2ii ~ 100 kHz) 
and Rydberg atoms whose transition frequencies can be a factor 10® 



28 



smaller 

The expression in eqn. ( |5.18[ ) is easily evaluated numerically. To properly eliminate the exponentially 
small but nonvanishing overlap of (^o(i') with the surface, we cut off the k and ^ integrations at suitably 
large values, as explained in detail in appendix [B] This procedure applies in exactly the same fashion to 



the integration in eqn. (6.91 below 



Figure [T] shows the energy shift of a rubidium atom in the harmonic trapping potential eqn. (5.4 1 
with i//27r = 1 kHz. At this frequency, the oscillator length is oq ~ 340 nm. As the overlap of the atom 
wavefunction with the surface should be negligible, we restrict the evaluation to a distance range d> 2 fim, 
making an error of the order Exp[-{d/aof] ~ 10"". The black lines in fig.[l] are for the case of a perfectly 
reflecting surface, with the reflection amplitudes R^ = 1 and R" — —1. The red lines involve a frequency- 
dependent reflection, as appropriate for a gold surface (described by the Drude model, see appendix [A| for 
details). The two terms of eqn. ( 5.18 1 are shown separately, the recoil correction (dashed lines) is multiplied 
by a factor of — cjeg/j^ to fit on the scale. The dashed horizontal line shows the asymptotic expression for 
the Casimir-Polder potential of an atom in front of a perfect mirror. 



AEcpid) 



cq(0) 3 
' eod" 8^ 



(5.22) 



with the static polarizability a(0). 



6 Ideal Bose gas in a surface trap 

Now we will use the results obtained in section [4. 2. 2| to calculate the atom-surface interaction for a gas of 
N noninteracting bosons. The atomic system (treated in the grand-canonical ensemble) is supposed to be 
in thermal equilibrium at the inverse temperature P, but its temperature is allowed to be different from 
the field temperature. 



6.1 Atomic correlation function 

As above in section[5] we will assume an isotopic harmonic trapping potential of the form eqn. ( 5.4 1. For the 
ideal Bose gas, and that enter the expression eqn. (4.151 for 'i/g are then simply the eigenfunctions 
and energies of a three-dimensional harmonic oscillator: 



-En = {rix +ny + n^) u, m = 0, 1, 2, 



(6.1) 



where the ground state of the trap has been set equal to the zero of energy. The critical temperature takes 
the value (25l[29] 

1/3 



(6.2) 
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d//im 

3 4 5 7 9 



< 
UJ 

<l 

I 



2.5 



2.25 



1.75 




Figure 1: Atom-surface interaction energy between a rubidium atom and a perfectly reflecting (black/upper 
solid) and gold surface (red/lower solid), field at zero temperature. Distance d — xc/uj^^g in units of the 
resonance wavelength, energy multiplied by and scaled by the natural linewidth h'jf.g — j/x'^^pCLig /(STreoc"^)- 



The atom is trapped in the harmonic potential eqn. (5.4) with a trap frequency 1^/27: — 1 kHz. Dashed 



lines: recoil correction multiplied by —ujeg/iy (see eqn. (5.18l). Resonance frequency ujeg = 27r3.85 x 10^'' Hz 



(isotropic polarizability); parameters of the Drude dielectric function for gold, eqn. (A. 5 1 



ujpT = 5 X 10^. 



= b.TAujeg and 
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2 



6/ium 

3 4 5 7 9 



2.7 



< 

I 



2.5 



2.3 




Figure 2: Atom-surface interaction between a trapped rubidium atom and a perfectly reflecting surface (field 
at zero temperature). Energy scaled by d"* as in fig. [l] The full curves correspond to a perfectly localized 
atom (black/bottom) and delocalized atoms with different trapping frequencies (y jlix — 3, 2, 1 kHz in blue, 
green, red from bottom to top). Top/pink curve: interaction energy per atom for a trapped ideal Bose gas 
at T = 0.2 Tc and v j^ix — 1 kHz (see eqn. ( |6.9[ )). Horizontal dashed line: asymptotic expression AEcp (see 
eqn. (5.22)) for the Casimir-Polder potential. Dashed curve: AEcp multiplied with the enhancement factor 
eqn. (|6.12| for !^/27r = 1 kHz. 
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For a given mean particle number A'' and inverse temperature /3, the (negative valued) chemical potential 
^i{N, (5) has to be determined from the relation 

iV(M,/3) = |dV{*t(r)*,(r)) , (6.3) 

where the brackets (...) denote a state of the atomic system at temperature p. 

With the same arguments as in section [53] above, the propagator for the excited atomic state is local 



(6.4) 



to a very good approximation. Neglecting the small recoil correction, we find from eqn. (4.181 



(*],(a;2)*s(xi))*,(2^2)*l(a;i) 



$n(r2)$n(ri) 

exp[/3(K - At)] - 1 

xe(t2 -ti)5(r2 -ri) . 



E 



exp[— i(t2 — t\){LOe 



En 



To the same precision, we can neglect the single particle energy .En compared to the atomic transition 
energy ujeg. This is even true for realistic atom temperatures: for a trapping frequency z//27r = 1 kHz and 
a mean number of TV = 10^ trapped particles, the mean thermal energy that sets the scale for the relevant 
En evaluates to T = (T/Tc) 2n 20.3 kHz much smaller than ujeg. We thus neglect En in the exponential in 
eqn. (6.41 and obtain 



e{t2 - h) 5{r2 - n) {*t (r2)*3(ri) )e 



-iuJeg(t2—tl) 



(6.5) 



The correlation function (\['g(r2)5'g(ri)) that enters in eqn. (6.31 and in eqn. (6.5l above reads 25 29 



{*t(^2)*,(^i)) = {v^aoy'J2\^'^"i^-^ 

1 , 



-2j73i'\-3/2 



(6.6) 



X exp 



^yr? + r?l^tanh[jM2]+lr^-rfl 



coth[j/3/./2]) 



where the vectors r'* = (x, z — d) account for the distance d between the surface and the center of the trap. 



6.2 Surface-induced energy shift 

With these approximations, the general expression eqn. ( |3.4[ ) gives a T-matrix 

duj Dai3{r,r, 



(TVIT^^'ITV) =McM/5 / dV(*t(r)*,(r)) 



(6.7) 



2n U} — uJeg + ie 

Had we kept the trap eigenenergy En, it would appear as a small shift of cjeg in the denominator. Using 
eqn. (4.3 1 and neglecting any thermal photons (see the remark below eqn. (4.3 1), we obtain 



{N\T'-^^\N) = ^^^lip- I d^r{*;(r)*3(r)) / d^G{r,r,t^) 



o^'ig+e 



(6.8) 



Performing the spatial integration and switching to the dimensionless variables of eqn. (5.181 yields 



(iv|r<^'liv) 



_L ^° j/3fi 

i=i 
, kdk 



tanh[j/3i^/2] 



-3/2 



_ I(Ji,x,rj+)Mc,p{k,^,R'',R') ^, 

1 + C 



(6.9) 



where I{K,x,ri+) is defined in eqn. (5.191, and the matrix Map in eqn. (5.9 1. The Lamb-Dicke parameter 
rj^ = a+ujeg/c now involves the temperature dependent width 

a+ = ao (tanh[j/3i//2])"^/^ > ao . (6.10) 
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To compare eqn. ( 6.9 1 with the resuh eqn. ( 5.18 1 for the single atom, we note that the constraint eqn. ( 6.3 1 
leads to 

oo 

e-"'f' ((1 - e-^'^-") tanh[j/3i//2])~^/'' = iV (6.11) 

and consider an interaction energy per atom, {A^IT'^' |A'')/A''. The terms with large j in the sum involve 
a width a+ equal to the zero-temperature value oq. These terms describe the condensate atoms in the 
trap ground state. The terms with small j have larger values of a+ and contribute to the energy shift as 
a broader trap would do. Indeed, for j = 1 and <^ 1, one gets the spatial width of a classical, thermal 
density distribution. 

This behaviour is shown in the numerical evaluation of eqn. ( 6.9 1 and eqn. (5.181 in fig. [2j for a perfectly 



reflecting surface. (More realistic materials can be described without further complications.) The atom- 
surface interaction per atom at T = 0.2 Tc is larger than for a single atom (at the same trap frequency 
h'/2n — 1 kHz), which is due to the larger spatial size of the thermally excited trap levels. 

At an atom-surface distance of d > 2/im, the interaction potential for the perfectly localized atom 



(calculated from eqn. (5.17l) is already deep in the retarded a;~^ regime. For an atom delocalized in the 
trap, the interaction potential becomes larger in magnitude because of the curvature of the Casimir-Polder 
interaction. Averaging a power law 1/z* over a narrow distribution {a <C d) centered at z = d, we get to 
leading order the enhancement factor 

(i>^l[l + 5(a/d)^ + ...] . (6.12) 



z 



The dashed black curve in fig.[2]shows the asymptotic expression for the Casimir-Polder potential eqn. ( 5.22 1 



multiplied with the above enhancement factor for a trapping frequency of v/2n = 1 kHz. The estimate 



eqn. (6.121 is seen to be in good agreement with our result from eqn. (5.181 (red line) 



7 Summary and outlook 

The starting point of our calculation was a second-quantized Hamiltonian that describes the interaction 
of a trapped system of A'^ atoms with the electromagnetic field. We have focused on two simple models 
for the atomic system: an interacting BEC described by A'^o atoms populating a single condensate wave 
function (described by the state |A''o)) and a noninteracting Bose gas at finite temperature, where the 
N particles populate the various single particle states of the trap (this state is denoted schematically by 
|A'^)). To calculate the interaction energy between the atoms and a plane surface, we made a perturbative 
expansion of the electromagnetic self-energy and worked out the T-matrix elements (7Vo|T'^^ | A*'o) and 
^jYlr'^' |A'') to second order in the atom-field coupling. The methods developed here are general enough 
to push the diagrammatic expansion to higher orders. The electric field propagator has been expressed 
in terms of retarded Green functions that permit to identify easily the contribution brought about by 
the surface. The characteristics of the surface material then enter through the scattering amplitudes for 
light, which allows for treating a wide range of materials. For the sake of simplicity, we considered the 
field to be at zero temperature as well, but thermal corrections can be included in a straightforward way 
by considering the temperature dependent term in eqn. (|4.3[|. Even non-equilibrium situations (bodies at 



different temperatures) can be covered by combining the techniques of fluctuation electrodynamics 30 
with the Keldysh formalism (see Ref. [3l] for an example). 

The expression found for (A^oIt'^' |iVo) in eqn. i5.2l describes the Casimir-Polder like interaction energy 



of a trapped Bose gas with the surface, for a general condensate wave function 0o(x). If the system is 
reduced to a perfectly localized single atom as treated in ^10, , our expression reproduces known results (see 



eqn. (5.17l). It also highlights that in full generality, the atom-surface interaction does not reduce to an 



integral over the density distribution of the atoms, due to the (virtual) propagation in the excited state. The 
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Bose gas-surface interaction energy shows an overall scaling with the atom number A^o (as can be expected 
at this order of perturbation theory), but even the interaction energy per atom still depends weakly on 
A'^o- We have identified for this dependence the following physical mechanisms, (i) The interaction energy 
involves a spatial average over the density profile whose width is larger for repulsive atom-atom interaction. 
This effect was already taken into account in the pioneering experiments of Ref. [5][7|. (ii) The atomic 
interactions (treated here as a contact potential) shift the optical transition frequency (see for example the 
experiments of Ref. [32]) and modifies the ground-state polarizability eqn. ( |5.13| |. (iii) The optical spectral 
line is recoil-broadened due to the kinetic energy of the atoms. This effect is very weak for typical traps 
and in the fully degenerate limit as the phase gradient of the condensate wave function vanishes. 

For the ideal Bose gas from (iV|T(^'|iV) in eqn. i6.9f, the Casimir-Polder interaction per particle does 
not depend on the atom number. We showed that the influence of a higher atom temperature on the 
atom-surface interaction is similar to that of a broadening of the trap potential. 

We plan to generalize the method presented in the present paper in two directions: on the BEC side of 
the problem, we want to include contributions from higher collective modes (condensate depletion, phase 
fluctuations, thermal density fraction) and revisit the problem of two atoms in front of a surface |33| . 
This setting has also been realized in a many-body version, by splitting a BEC in two spatially separated 
modes 34,35 . On the cavity QED side, higher orders of atom-photon interactions will be considered where 
intensity fluctuations of the quantum field near the surface appear [36|l37|. 



A Retarded Green function for the electric field in the pres- 
ence of an interface 

The reflected part of the retarded Green function in the presence of an interface (as presented in [oj 
eqn. (3.4)], see also [38[ sec. 2] for an overview) reads 

/-ift I \ f ^ ^ n /^ \ ifc,(2i +22)+ik.(xi — xq) /a i\ 

Ga/3(ri,r2,aj) = ^ 2^600^ Rc0(k,uj) e (A.l) 

with kz = \/ oj^/c^ — Here, the two-dimensional vectors x and k denote the position and momentum 
vectors parallel to the surface, respectively. Henceforward in this appendix, we use units with c = eo = 1. 
The matrix Rap (k, a;) is deflned as 

i?Q/3(k, tj) = (ss)aa R" -f (po+po-)a/3 R^ (A. 2) 



The functions R^ and W in eqn. (A.2l are the Fresnel reflection coefficients for s- and p-polarized light, 
which can be modeled to realize different surface materials. For the case of a perfectly reflecting surface, 



R" = —\ and RF = \, while in general the reflection coefficients are frequency dependent (see [9j[l0 24 ): 
Considering an interface between vacuum (eo = 1) and a material with a local and isotropic dielectric 
function R" and R^ axe given by 

fe, -f (a;2^(a;)-fc2)i/2 ' 

^ ekz-{u;'e{uj)-kY^' 
efc, + (a;2e(a;)-fc2)i/2 ' 

In section [STI} we use the Drude model for a metal surface, with 

2 

e{Lo) = 1 - , , , (A.5) 
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where cjp is the plasma frequency and r the collision time. Finally, the dyadic elements (ss)q,^ and 
(po+Po-)c«/3 in eqn. (A.2I involve the normalized polarization vectors 



Po± = 



— k X z 

z =F fc^k 



(A.6) 
(A.7) 



B Approximating the error function integral 

In the integrands of eqn. | |5.18[ ) and eqn. ( |6.9[ ), we encounter the expression 

1 X 

I{k, x,rj) = - exp[— 27t2: + k r; ] (1 + erf [ Kr)]) (B.l) 

2 r? 

where k is integrated from zero to infinity, rj = aoUeg/c is fixed by the atomic transition frequency and 
mass and the trap geometry and the positive distance x varies such that a; > r;^ is always fulfilled. 

Noting that the argument of the error function changes sign a.t k — x/rj, we can approximate the error 
function for large values of k (see [39| eqn. (8.254)]) to obtain 

rr- \ exp[— a;^/j7^1 r - ^ /-n 

I{k,x,t^) ^ ' ' ' , forK»— B.2 

which is exponentially small in the quantity {xjr])^ . In numerical integrations, we will thus cut off the 
ci«:-integration at k = x/rf, omitting terms of order C'(exp[— (a;/?;)^]) in the integrand. The neglected 
quantities are small: for a rubidium atom at T = trapped in a v/2n = 1 kHz trap at an atom-surface 
distance of a; = ducg/c = 15, we have {x/rff ~ 30. Conceptionally, the high momentum cut-off is necessary 
as the atomic probability density |(^o(i')P we adopt here is not zero at the surface, but only exponentially 



small, namely of the same order as the terms neglected in eqn. (B.2 1. 
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